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disorder are investigated. We employ the statistical dynamical mean- field theory, which incorporates 
both, fluctuations due to disorder and local correlations due to interaction, to solve the Anderson- 
Hubbard model. Localization due to disorder is studied by means of the probability distribution 
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diagram consisting of disordered correlated metal, Anderson-Mott insulator, and band insulator. 
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I. INTRODUCTION 

The interplay between disorder and interaction in cor- 
related electron materials still remains far from com- 
plete understanding^ in spite of major progress achieved 
in disordered but non-interacting electron systems 
On the theoretical side, the main obstacle is the non- 
perturbative character of the most interesting phenom- 
ena, e.g., alloy-band splitting, Anderson localization, or 
the Mott-Hubbard metal-insulator transition.^ On the 
experimental side, the main problem is that both the 
amount of disorder and the strength of the interac- 
tion are not well controlled and not easily tunable in 
real condensed matter systems. The idea of a quan- 
tum simulator, 5 where a complicated quantum many- 
body system is simulated by another quantum but per- 
haps simpler system working as a quantum computer, is 
very attractive. Regarding the interplay between disor- 
der and interaction, experiments with ultracold fermionic 
or bosonic atoms in optical lattices^H' are very promis- 
ing steps toward creating such a quantum simulator and, 
therefore, are capable of shedding light into this unsettled 
problem. 

Disorder in ultracold gases can be simulated in dif- 
ferent ways: (i) by using an optical speckle laser Q^EHl 
(ii) by superimposing two laser beams with incommensu- 
rate frequencies,^ or (iii) by loading two atomi c spe cies, 
where only one is mobile, into an optical lattice! 15 * 16 ! The 
latter simulates a binary-alloy distribution of the on-site 
energies. Effects of interactions in optical lattices are 
controlled by tuning the on-site potential depths and/or 
the magnetic field around a Feshbach resonance.^ Re- 
cently developed momentum resolved radio frequency 
(rf) spectroscopy,^ on the other hand, is an adequate 
probing technique of correlated and disordered systems. 
This technique is similar to the well-known angle resolved 
photoemission spectroscopy experiments j 19 l 2Q l which were 
performed for alloys in solid-state physics. rf spec- 
troscopy probes the spectral function and, thereby, the 



single-particle Green's function of the many-body sys- 
tem. 

Strongly correlated fermions in three dimensions are 
successfully described within the dynamical mean-field 
theory (DMFT)P^ This mean-field theory is fully 
non-perturbative and in combination with density func- 
tional theory it is capable of describing properties of 
real solid-state systems.^ An extension of DMFT, which 
includes disorder effects, was performed both, in anal- 
ogy to the w ell-known coherent potential approximation 
(CPA)p^^and within a fully stochastic appro ach to in- 
corporate effects of Anderson localization ! 27 * 28 ^ Since the 
latter approach is computationally very expensive if one 
treats correlation effects on a rigorous level and keeps 
sufficiently large ensembles of disorder realizations, the 
typical medium theory (TMT-) DMFT was developed. 1 ^ 1 
Here the geometrically averaged local density of states 
(LDOS) is used as an order parameter for Anderson lo- 
calization. TMT-DMFT was successfully applied to the 
non- interacting 2 ^ and to the interacting^"^ electron sys- 
tems with disorder. 

By construction, TMT-DMFT is only capable of de- 
scribing effects of strong localization due to disorder, i.e., 
effects caused by fluctuations of the wave- function ampli- 
tudes. Since TMT-DMFT determines the typical LDOS, 
i.e. the most probable value of the LDOS, all non-local 
phase interference effects are missed. To improve the 
theory such that weak- localization effects in the many- 
particle wave function are kept, one should combine a 
fully stochastic approach with DMFTP^In this way one 
can readopt the original point of view of Anderson 2 and 
use the full probability distribution function (PDF) of 
the LDOS as an order parameter for the Anderson tran- 
sition within DMFT. Recently, also the periodic Ander- 
son model with disorder was investigated by means of the 
statistical DMFTP^MD Therein, a novel electronic Grif- 
fith's phase, characterized by non-Fermi-liq uid behavior, 
was established as a precursor of a disorder-driven metal- 
insulator transition. For non-interacting disordered sys- 
tems such a stochastic theory, named local distribution 




(LD) approach, was effectively used in an analytical 
approach^ and recently implemented numerically.^ 

The aims of this paper are to apply the statistical 
£)]VIp r l[2l|38] ^ interacting and disordered fermions and 
to extend the method to a level not reached so far. 
The theory is applied to correlated fermions on a lat- 
tice with binary- alloy type of disorder. This problem has 
recently been addressed within DMFT combined with 
CPA to deal with disorder.^ It was shown, in particu- 
lar, that new types of alloy-Mott or alloy-charge trans- 
fer insulators can appear and that the Mott-Hubbard 
metal-insulator transition can occur at non-integer parti- 
cle densities! 3 -^ Here we revisit this model and show that 
Anderson localization significantly extends the picture. 
In order to make the statistical DMFT method com- 
putationally feasible the DMFT part is solved approxi- 
mately within modified perturbation theory (MPT) S^IMJ 
which (in c ontrast to, e.g. the slave boson mean- field 
th eor -^23]42|^ provides a reliable interpolation scheme be- 
tween the weakly and strongly interacting regimes.^ 
However, a quantitative analysis of the above mentioned 
Griffiths phase is limited by MPT, as this impurity solver 
does not reproduce the exponentially small low-energy 
scale for strong interactions. 

The paper is structured as follows: In Sec. [H] we 
introduce and motivate the underlying physical model. 
The LD method for non-interacting disordered systems 
is reviewed in Sec. IIII Al and extended to the statistical 
DMFT for interacting systems in Sec. |IIIB Our main 
results are discussed in Sec. [IV| Finally, the connection 
to experiments in optical lattices is discussed in Sec. [V| 



II. ANDERSON-HUBBARD MODEL WITH 
BINARY-ALLOY DISORDER 

Electrons or cold fermionic atoms, such as 6 Li or 40 K, 
in disordered lattices are well described by the Anderson- 
Hubbard Hamiltonian 
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where c\ a (ci a ) denotes creation (annihilation) operators 
at a lattice site i with spin a = ±1/2. The fermionic 



number operator is given by 



Ci a . The hopping 



amplitude between sites i and j is denoted by Uj, the 
interaction amplitude is represented by U, and the chem- 
ical potential is given by /i. In the following we consider 
fermions on a Bet he lattice^ with connectivity which 
is related to the coordination number Z via K = Z — 1, 
where the hopping amplitude is only nonzero = t for 
nearest neighbors i and j. We also set energy units such 
that the band-width Wo = ^t\fK — 1 hereafter. The lo- 
cal disorder is given by random on-site energies e^, which 
are drawn from a probability distribution function p e (ci). 

In this paper we consider the case of a binary-alloy 
Anderson-Hubbard model, in which the PDF of the on- 



FIG. 1: (Color online) Illustration of a realization of binary 
disorder in optical lattices. Two atomic species (indicated as 
light green and dark blue spheres) are loaded into an optical 
lattice. The hopping amplitude of one species (dark blue) is 
suppressed and therefore these atoms are immobile. Due to 
the interatomic interaction the second species experiences a 
binary disordered lattice potential depending on the presence 
of an atom of the immobile species on the same lattice site, 
i.e. if there is a dark blue atom present the on-site energy is 
€i — +A/2 otherwise 6i — —A/2. 



site energies is given by the bimodal function 



p e (ei) = xS(ei + + (1 - x)5(e i - ^), (2) 



x are the fractions of lattice sites with 

A 



where x and 1 

energies = — y and = respectively, and A de- 
scribes the on-site energy splitting. In general, A and x 
are independent parameters. However, the cases x = 
or 1 correspond to non-disordered systems with on-site 
energy shift ±A/2. Therefore, a natural parameter for 
measuring the disorder strength in binary alloy systems 

is S EE X(l ~ X)AP 

A very important difference between binary-alloy dis- 
order and disorder types with continuous probability dis- 
tributions is that in the former case in a non-interacting 
system and in arb itrary lattices the Bloch band is split 
if A > Wq! 39 ' 45 ' 46 1 In this limit two alloy subbands are 
formed and the system is a band insulator if v = 2x or 
v = 2, where v is number of fermions per site, or a metal 
otherwise. In the presence of interaction a Mott insula- 
tor at fractional particle filling ^ = xor^ = l + xis 
allowed! 39 ^ 47 ! Here we investigate how Anderson localiza- 
tion modifies these predictions. 

In systems of cold atoms in optical lattices the binary- 
alloy disorder is prepared by adding an additional species 
of atoms, which are immobile but interact with the mo- 
bile components. First experimental attempts in this 
direction have been performed! 15 * 16 * However, in such a 
system one must take care that the immobile atom po- 
sitions are random but not fluctuating in time, i.e. the 
created disorder must be quenched. 48 Such a situation is 
schematically presented in Fig. [I] 



III. METHOD 

In this section we introduce the notation used and de- 
scribe the statistical DMFT for non-interacting and for 



Stochastic Bath 



interacting systems in Sees. Ill A and IIIB[ respectively. 



A. Local distribution approach 

The local distribution approach is a self-consistent 
computational scheme for determining the probabil- 
ity distribution function of the local single-particle 
Green's functions, i.e. p[Ga a (uj)\. Here Gij a (uj) is the 
Fourier transformation of the retarded Green's function 
G ija (t) = -iO(t)([ci a (t),c] a (0)]+), where 0(t) is a Heavi- 
side function and [..,..]+ denotes anticommutator brack- 
ets. In the following we consider only paramagnetic solu- 
tions of the Anderson-Hubbard model and therefore the 
spin index a is omitted. 

In the absence of interactions the renormalized pertur- 
bation theor}J^ shows that the local Green's function can 
always be expressed as 



Gu(u) 



(3) 



where the hybridization function Ti(uj) describes all ef- 
fects of the coupling of site i with other nearest neighbor 
lattice sites. The chemical potential is given by p. For 
numerical reasons we also introduced the broadening fac- 
tor r] > 0. In order to study localization effects, the limit 
r] — >• has to be performed. 

The hybridization function Ti(uj) can be expressed by 
an infinite, renormalized series of the form 



K 



(4) 



where G^} (uj) is the cavity Green's function of the system 
when the site i is removed. 

On the Caley tree (Bethe lattice) this series can be ex- 
actly truncated after the first term and the hybridization 
function is exactly given 



Ti(u) = t 



3 = 1 



(5) 
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FIG. 2: (Color online) Illustration of the statistical dynami- 
cal mean-field theory applied in this work. The many-body 
problem with disorder is mapped onto an ensemble of single 
impurities, which are coupled to an ensemble of stochastic 
Green's functions, which is determined self-consistently. G n 
represents the nth sample from the ensemble of Green's func- 
tions. 



single-particle Green's function Gu(u) is calculated using 
Eq. ([3]). (iv) Having calculated all new Gu(u) a new PDF 
p[Gii(uj)] is obtained and we return to step (i). The al- 
gorithm is repeated until self-consistency for p [Gu(oj)] is 
achieved. We note that this method incorporates spatial 
fluctuations, i.e. quantum interference effects, caused by 
the disorder. Schematically the computational procedure 
is presented in Fig. [2j 

The relevant physical observable is the LDOS Pi(oj) = 
— ^lm(Ga(uj)) : which is a random quantity in disordered 
systems. The corresponding distribution p[pi(uj)] is ob- 
tained by counting all values of the LDOS for each fre- 
quency and constructing a histogram.^ From this proba- 
bility distribution we can then determine the expectation 
value, i.e. the arithmetically averaged LDOS 



(pM)arith = (PiM)dis 



(6) 



and the typical value, which we approximate by the geo- 
metrical average 



(pM)geom = exp(lnp i (o;)) dis 



(7) 



with K as the coordination number of the lattice. Here 
we use the fact, that the geometry of the Bethe lattice 
does not change when site i is removed, which allows us 
to determine the cavity Green's function G^ in analogy 
to Gu. 

In practice, given an initial PDF p [Gu(uo)) the compu- 
tational scheme is the following: (i) For each ensemble 
member we draw a random on-site energy out of the 
PDF p e (ei) given in Eq. (ii) The hybridization func- 
tion Ti(uj) is determined via Eq. ([5|, in which the near- 
est neighbor cavity Green's functions G^ (lj) are ran- 
domly sampled from the PDF p[Ga(oj)]. (hi) The local will also be useful to characterize the disordered system. 



where (F) dis = J °° dxF(x)p[x] is the average over differ- 
ent disorder realizations of the corresponding quantity F. 
In the following, the cumulative probability distributions 



p(u>) 

>[p(uj)} = J p[pf(uj)]dpf(uj) 



(8) 
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B. Interacting systems 

In the presence of interactions, Eq. Q is no longer 
true. It should be replaced by a Dyson-like equation, 
which relates the inverse of GV/(cj) -1 with the self-energy 
functions E^cj). Within the statistical DMFT the full 
self-energy E^ (u) is approximated by a self-energy diag- 
onal in the lattice indices, i.e. E^ 



Within 



this approximation for interacting systems, Eq. Q 
modified to 



1 



(9) 



This approximation of the self-energy becomes exact 
in infinite dimensions, as was shown by Metzner and 
VollhardtP^, and was u sed as a starting point for devel- 
oping the DMFTPEH 

Explicitly, within the statistical DMFT Hubbard 
model is mapped onto an ensemble of Anderson sin- 
gle impurity models (as schematically shown in Fig. 2]). 
We now repeat the previously described LD algorithm 
with the additional description on how to determine the 
self-energy of the interacting system. 

In a fully interacting problem, different frequencies uo 
do not decouple in the self-consistency relations. This is 
in contrast to the non-interacting case, where the self- 
consistency equations are solved for each frequency sepa- 
rately. In the interacting case, we are therefore restricted 
to ensembles typically of the order 10 3 samples. In or- 
der to reduce computation time, we also need to use a 
fast method (a so-called impurity solver) for determining 
the self-energy E(gj). Here, we use the iterative pertur- 
bation theory (IPT) ; H5l 4i | 53 | w j 1 i c j 1 properly reproduces 
the non-interacting and atomic limits and was shown to 
qualitatively describe the Mott-Hubbard metal- insulator 
transition at a critical interaction Within the IPT 
the self-energy is calculated in second order in U in the 
non-renormalized perturbation expansion. 

The original formulation of the IPT was restricted to 
the half-filled case. Later the method was extended to 
densities a way f rom half-filling; this is commonly referred 
to as MPTp^The self-energy within MPT is given 



E(cj) = Un 
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where 



and 



n(l — n) 
n(°)(l-n(°)) 



B-B^> -fjL + fi + U(l -2n) 
[/ 2 n(°)(l-n(°)) 



(10) 



(11) 



(12) 



are additional coefficients in the interpolative formula. 
Here, denotes the filling obtained by using the 

Hartree-Fock solution 



1. 



Jl — e — Un — T(uj) + in ) 



(13) 



The parameter Jl and the higher order correlation func- 
tion and B have to be fixed in such a way that 
the correct first three moments of the spectral density 
are guaranteed. In perturbation theory the second-order 
contribution to the self-energy is given by 



C/2 

i 



dt exp(iwi) (p- {t)p+ {t)p+ (t) 



+p+(-t)p-(-t)p-(-t)) , (14) 
where the Laplace transformed density of states is 



oo 

p ± (t) = J dwexp(-iujt)p ( ' 0) (±uj) . 



(15) 



According to Potthoff et alW^ there are three ap- 
proaches for fixing Jl. The first is to require \i — Jl. 
Second, one imposes the Friedel sum rule to ensure the 
low energy Fermi liquid behavior as done by Kajueter and 
KotliarP^The last possibility requires that = n. All 
three possibilities do not affect the validity of the MPT 
in the weakly interacting limit as all methods guaran- 
tee that jl — » fiu=o as U — » 0. 40 Furthermore, all three 
approaches have been compared carefully and checked 
against exact diagonalization (ED) calculations.^! In con- 
clusion, the second and the third approaches show very 
good agreement whereas the first one differs considerably 
from ED results. In this work we choose the third possi- 
bility. 

The higher order correlation function B^ is expressed 
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xG (0) (cj + ^). 
The correlation function B is given by 



Im / dujY(uj + in) 



(16) 
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dujT(uj + in) 



x (^£(cj + irj) - 1)G(cj + in) , 



(17) 



which can be solved self-consistently. 

Studying correlated and disordered lattice fermions 
within statistical DMFT, we found that it is useful to 
compare the resulting spectra l functions to those deter- 
mined within the CPA ESEHEI Within CPA the hybridiza- 
tion r(cj) is given by 

r(w) =t 2 K[xG tz= _ A/2 (co) + (l-x)G tz=+A/2 (co)}. (18) 

Note however, that the CPA is not able to describe An- 
derson localizationP^^ 
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FIG. 3: (Color online) Natural logarithm of the PDFs of 
the non-interacting system p[p] plotted color coded for fixed 
broadening r\ = 10 -3 and several disorder parameters A: (a) 
A = 0.0, (b) A = 1.0, (c) A = 2.0, (d) A = 3.0. Parameters 
are K = 6, v = 0.2, x 0.2. 
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FIG. 4: (Color online) Geometrically (dashed red line) and 
arithmetically averaged (solid black line) spectral functions 
of the non-interacting system for fixed broadening 77 = 10 -3 
and several disorder parameters (A = 0.0,1.0,2.0,3.0,4.0). 
Parameters are K = 6, v — 0.2, x — 0.2. 



IV. RESULTS 

In this section we present our results obtained by sta- 
tistical DMFT concerning Anderson and Mott- Hubbard 
transitions in correlated fermionic systems with binary- 
alloy disorder at zero temperature. In particular, we in- 
vestigate how the predictions from Ref. [39] are extended 
when the Anderson localization is present. 

In the following we set the impurity concentration x 
and the total particle density v to be equal, i.e. x = z/, 
by adjusting the chemical potential during the iterative 
solution of DMFT equations. This choice enables us to 



study Mott-Hubbard metal-insulator transition at non- 
integer particle densities. For practical calculations, we 
choose the impurity concentration and the particle den- 
sity equal to x = v = 0.2. Furthermore, we set the 
coordination number K = 6 shortly above the classical 
percolation threshold x p = 1/ifp^i.e., extended states 
can exist within both upper and lower alloy bands when 
they are split due to disorder. 



A. Detecting Anderson transition in the 
non-interacting case 

We first discuss how to detect localization effects and 
how to distinguish between extended and localized states 
in the non-interacting limit. Figure [3] shows the PDFs 
on a logarithmic scale for different increasing values of 
the disorder parameter A. A band splitting into an up- 
per and a lower alloy band occurs with increasing dis- 
order. This can also be seen plotting the arithmetically 
and geometrically averaged density of states as shown in 
Fig. [4j Looking at Fig. |4j we observe a vanishing geo- 
metrically averaged LDOS in the minority band. This 
corresponds to disappearing extended states and is used 
within TMT-DMFT to identify the Anderson transition. 
As mentioned before, in this paper we use a more power- 
ful and general approach to detect Anderson localization. 

Extended states are characterized by a branch cut on 
the real axis of the local Green's function, whereas lo- 
calized states are characterized by a dense distribution 
of poles in the thermodynamic limit This fact can be 
used to detect if states are localized or extended by inves- 
tigating the behavior of the PDFs of the LDOS p[pi(uj)] 
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FIG. 5: (Color online) Natural logarithm of the PDFs p[p] 
of the non-interacting system plotted color coded for fixed 
disorder parameter A = 3.0 and several broadenings 77: (a) 
ry = 10~ 2 , (b) j] = 10" 3 , (c) j] = 10" 4 , (d) 77 = 10" 5 . Param- 
eters are K — 6, v — 0.2, x — 0.2. 
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FIG. 6: (Color online) Comparison of color coded natural 
logarithm of PDFs p[p(u>)] of the minority band of the non- 
interacting system for two different disorder parameters A = 
1.0 (plots on the left side (a)-(d)), A = 4.0 (plots on the right 
side (e)-(h)) and for several broadenings 77: (a,e) 77 = 10 -3 , 
(b,f) 77 = 10" 4 , (c,g) 77 = 10" 5 , (d,h) 77 = 10" 6 . Parameters 
are K = 6, v = 0.2, x 0.2. 



shown in Fig. [3] when the broadening 77 tends to 0P^ 
Namely, the PDF of the LDOS for extended states satu- 
rates at a finite value for 77 — >> 0, whereas the PDF of the 
LDOS for localized states decreases to zero for 77 — > 0. As 
an example, Fig. [5] shows the behavior of the PDF when 
decreasing the broadening from 77 = 10 -2 to 77 = 10 -5 
for a selected value A = 3.0. A change is seen for states 
in the lower alloy band, whereas the PDFs of the up- 
per alloy band remain almost unchanged in this regime 
of 77. The probability distributions of the LDOS of the 
lower alloy band are presented in detail for A = 1.0 and 
A = 4.0 in Fig. [6j It is clearly visible that the PDFs for 
small A = 1.0, corresponding to the lower alloy band, 
become 77-independent for 77 — » 0. On the contrary, at 
large A = 4.0 the PDFs strongly depend on 77. 

In addition to Anderson localization effects we also ob- 
serve that the spectrum is fragmented (cf. Figs. [4] and 
[6J, due to the presence of states with different physical 
properties. These states differ in the behavior of the PDF 
of the LDOS for 77 -> (cf. Fig. [6} and are identified 
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FIG. 7: (Color online) Behavior of cumulative PDFs P[p(u)] 
of the non-interacting system with decreasing broadening 77 
(a) for an extended state at A = 1.0 and oj = 0.0, (b) for an 
Anderson localized state at A = 4.0 and oj = 0.0 and (c) for 
an "anomalous" localized state at A = 4.0 and 00 = —0.01. 
Parameters are K — 6, v — 0.2, x — 0.2. 



either as cluster resonance d 61 * 62 * or as "anomalous" lo- 
calized states!^ The resonance states are similar bound 
eigenstates but with a finite life-time. They appear be- 
cause of special geometrical configurations of the impu- 
rity atoms. The "anomalous" localized states are in fact 
extended states over the whole lattice but they are insu- 
lating and do not contribute to the dc conductivity.^ On 
a bipartite lattice these states have small wave function 
amplitudes on one sublattice and large amplitudes on the 
other sublattice. The typical 77 behaviors of the PDFs for 
given frequencies are shown in Fig. [7[ Panel (a) shows 
the behavior of an extended state, panel (b) presents the 
behavior of an Anderson localized state, and panel (c) 
shows the behavior of an "anomalous" localized state 
with its typical bimodal structure.^ 



B. Anderson and Mott transitions in the 
interacting case 

In the interacting limit we restrict our investigation of 
the 77-dependence to the lower limit 77 = 10 -5 , as we use 
small ensembles due to computational limitations. We 
also note that the MPT requires a small finite broadening 
in any case. 

We first comment on the denning properties of the dif- 
ferent phases arising. The paramagnetic metal is charac- 
terized by a non- vanishing arithmetically averaged LDOS 
at the Fermi level (p(ct/)) arith . The paramagnetic metal is 
gapless and hence compressible. Since it is also different 
from the Anderson- localized phase, the geometrically av- 
eraged LDOS at the Fermi level is finite as well. The 
Mott insulator possesses an excitation gap which is of 
the order of the interaction strength, and therefore this 
phase is characterized by a vanishing arithmetically av- 
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FIG. 8: Phase diagram for the interacting and disordered 
system in A — ?7-plane, showing of Anderson-Mott insulator, 
paramagnetic metal and band insulator. The solid line cor- 
responds to the transition between insulating phases and the 
metal, the dotted line corresponds to a vanishing arithmetic 
average of the LDOS at the Fermi level and the dashed region 
denotes the crossover between Anderson-Mott insulator and 
band insulator. Parameters are K = 6, v — 0.2, x = 0.2. 
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FIG. 9: (Color online) Arithmetic average of the LDOS for 
increasing interaction strength U (U = 0.0, 0.5, 1.1, 2.0, 2.5) at 
fixed disorder parameter A = 2.0. The results from statistical 
DMFT calculations (solid black line) are compared to CPA 
results (dashed red line). The inset shows the arithmetic mean 
of the LDOS at the Fermi level with increasing interaction 
strength. Parameters are if = 6, zy = 0.2, x = 0.2. 



eraged LDOS at the Fermi level (p(uj = 0)) arith . 

In the presence of both interaction and disorder an 
Anderson insulator with localized one-particle wave func- 
tions is not well-defined anymore due to many-body ef- 
fects. Therefore, we refer to an Anderson-Mott insulator 
phase if the PDF of the LDOS tends to zero at the Fermi 
edge uj = when r] — >• 0. 

We also need to distinguish between an Anderson-Mott 
insulator and a band insulator. The band insulator is 
characterized by (p(oo = 0)) arith = 0, but in contrast to 
the Mott insulator the excitation gap is determined by 
the energy distance between the upper edge of the occu- 
pied band and the upper alloy band, which in this case 
is proportional to A. 

The phase diagram presented in Fig. [8] is the main re- 
sult of our paper. We find a metallic phase which turns 
into a Mott insulator at small A due to alloy band split- 
ting and the mechanism described earlier in Refs. 1471 and 
l39l However, as the current results prove, this type of 
Mott-Hubbard metal-insulator transition at non-integer 
particle densities is also possible, if Anderson localiza- 
tion effects are taken into account. In the limit of large 
disorder parameter A the metallic phase as well as the 
Mott-Hubbard transition is terminated by Anderson lo- 
calization. The states in the upper part of the phase 
diagram in Fig. [8] are localized due to strong disorder. 

Spectra corresponding to the Mott-Hubbard transition 
are displayed in Fig. [9j where the arithmetically aver- 
aged LDOS obtained within statistical DMFT is com- 
pared to that obtained within a CPA type treatment of 
disorder. With increasing interaction U at fixed A three 
peaks emerge because of the Mott-Hubbard and band 
splitting transitions. Moreover, we observe additional 



spikes in the LDOS similar to those observed for the non- 
interacting system.^ These spikes are not reproduced by 
a CPA treatment of disorder and we conclude that they 
are due to local interference effects on clusters of impurity 
atoms. In the inset of Fig. [9j the arithmetic average of 
the LDOS at the Fermi level is presented as a function of 
U. The Mott-Hubbard transition appears to take place 
at U = 1.1. However, this is not a true transition point as 
it corresponds to the regime within the Anderson-Mott 
insulator where all states are already localized, cf. Fig. [8] 
We also see in Fig. [9] that by further increasing the in- 
teraction to U = 2.0 the upper alloy band and the up- 




FIG. 10: Arithmetic average of the LDOS with increasing 
disorder parameter A (A = 0.25,0.75,1.25,2.0,2.5) for fixed 
interaction strength U — 0.5 and broadening 77 = 10 -3 . Pa- 
rameters are K — 6, v — 0.2, x — 0.2. 
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electron-photon interaction <£. If the investigated system 
is a disordered alloy, the ensemble average of the pho- 
tocurrent has to be calculated.^ In Ref . [671 it was shown 
that the ensemble averaged photocurrent is essentially 
given by the Bloch spectral function which is defined ad^ 



A(cj,k) = Im(k|Tr(G(u;))|k) , 

7T 



(20) 
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FIG. 11: (Color online) Comparison of the color coded natural 
logarithm of PDFs p[p(u)] for disorder parameter A = 2.0 and 
several interaction strengths U: (a) U = 0.0, (b) U — 0.25, 
(c) U = 0.5, (d) U = 0.75. Parameters are K = 6, v = 0.2, 
x = 0.2. 



per Hubbard band are merging. This corresponds to the 
crossover regime between alloy Anderson-Mott insulator 
and alloy-charge band insulator indicated by the dashed 
area in the phase diagram in Fig. U 39 An additional ef- 
fect is observed in Fig. [9j namely, with increasing U the 
position of the upper alloy band is shifted with respect 
to zero on the energy scale, cf. Ref. |64j This shift of 
the upper alloy band resembles a situation seen in the 
exactly solvable atomic limit 

For comparison we present the LDOS at a selected U 

~ As 



value for different disorder parameters A in Fig. 10 



expected, we observe a band splitting with increasing A 
and the formation of an energy gap proportional to A 
between the lower and the upper alloy bands. Note that 
additional peaks appear in the lower band when the dis- 
order parameter A is increased. These peaks do not occur 
in a CPA treatment of disorder. 



Finally, in Fig. 11 we show the evolution of the PDFs 
across the Mott-Hubbard transitions. The onset of a 
three peak structure is seen as well as sharp resonances in 
the LDOS being broadened and washed out by increasing 
the interaction strength U. 



MEASURING THE DENSITY OF STATES IN 
SOLIDS AND ULTRACOLD FERMIONS 



In angle-resolved photoemission spectroscopy 
(ARPES) the photocurrent is given by 6 - 6 -! 

I ph (uj) = -hm J J drdr'^ 2 (r,o; + ft)<I>(r) 

xG(r,r / ,cj)<I>*(r / )^2(r / ,^ + ^), (19) 

with G as the single-particle Green's function of energy 
u;, the LEED state ^2, the photon frequency ^, and the 



where (...) denotes the ensemble average. Here, the mo- 
mentum k is a good quantum number as the ensem- 
ble averaging restores translational invariance. Hence, 
in ARPES applied to alloys the Bloch spectral function 
is measured and therefore the ensemble averaged single- 
particle Green's function, which allows a comparison of 
the experimentally measured density of states to theoret- 
ical calculations (cf. Ref. l69l and references therein). 

Regarding experiments with cold atoms, a very promis- 
ing probing technique, the momentum-resolved radio 
frequency spectroscopy,^ has been developed. In the 
experiment!^ a radio frequency field of frequency Q rf was 
applied to a two-component mixture of fermionic 40 K 
atoms in hyperfine states |1) and |2) in order to excite 
the atoms of hyperfine state |2) to hyperfine state |3). 
The trap is then turned off and by counting the number 
of atoms in state |3), A/3, the dispersion ek is obtained.^ 
The rf current, defined by I = (A3) is given b}E^ 



j(k,H = M p(k)W)/HL 



(21) 



for homogeneous systems. Tk denotes the transition ma- 
trix, 5v is the rf detuning, and £k is equal to k 2 /2m — /i, 
assuming that hyperfine state |3) is not occupied. In 
comparison, for homogeneous systems the photocurrent 
of Eq.([l9|) reduces tcP 



J ph (k,o;)=M(k,fi)p(k,a;)/(a;), 



(22) 



which, compared to Eq.(21), reveals the analogy between 



ARPES and momentum-resolved rf spectroscopy for ho- 
mogeneous systems. In inhomogeneous systems, such 
as trapped and/or disordered systems, final state effects 
have to be taken into account, which can be described 
by the density functional theory with the local density 
approximation.^ In conclusion, the arithmetically aver- 
aged LDOS calculated in this work can in principle be 
compared to spectra resulting from momentum-resolved 
rf spectroscopy applied to fermions in an optical lat- 
tice. In order to realize such a comparison, the analog of 
Eq. (19) needs to be calculated for the RF current and 



the ensemble average has to be carried out accounting 
for final state effects. 65 



VI. SUMMARY 

We have investigated the binary alloy disordered Hub- 
bard model within statistical DMFT, using MPT as an 
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impurity solver. This method treats disorder and inter- 
action on equal footing and in a non-perturbative way. 
The scheme reduces to the local distribution approach for 
non-interacting systems and to a standard DMFT-MPT 
scheme in the pure case. Applying the statistical DMFT, 
we were able to compute the full probability distribution 
function of the local density of states, and therefore, lo- 
calization effects have been studied in a more rigorous 
way and in more detail than in a typical medium the- 
ory combined with DMFT. As a result, the paramagnetic 
ground state phase diagram was obtained. It consists of a 
disordered metallic phase, an Anderson-Mott insulator, 
and a band insulator. For non-integer particle density 
n = x a Mott-Hubbard metal-insulator transition is ob- 
tained even when Anderson localization effects are taken 
into account. 

In future work, the method will be generalized to finite 
temperatures. It can also be extended to the hypercu- 



bic lattice. From the methodical point of view, it would 
be desirable to use impurity solvers that are superior or 
complementary to MPT. 



Acknowledgments 

We acknowledge many useful discussions with D. Voll- 
hardt and U. Bissbort. This work was supported by the 
Deutsche Forschungsgemeinschaft (DFG) via Forscher- 
gruppe No. FOR 801. Computations were performed at 
the Center for Scientific Computing of the Goethe Uni- 
versity Frankfurt. K.B. also acknowledges support by the 
Sonderforschungsbereich No. 484 of the DFG and Grant 
No. N202026 32/0705 of the Polish Ministry of Science 
and Education. 



1 N.F. Mott, Proc. Phys. Soc. A 62, 416 (1949); Metal- 
Insulator Transitions, 2nd edn. (Taylor and Francis, Lon- 
don 1990). 

2 P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

3 B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 
(1993). 

4 M. Imada, A. Fujimori and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

5 R.P Feynman, Int. J. Theor. Phys. 21, 467 (1982). 

6 D. Jaksch, C. Bruder, J.I. Cirac, C.W. Gardiner and P. 
Zoller, Phys. Rev. Lett. 81, 3108 (1998). 

7 M. Greiner, O. Mandel, T. Esslinger, T.W. Hansen and I. 
Bloch, Nature 415, 39 (2002). 

8 W. Hofstetter, J.I. Cirac, P. Zoller, E. Dernier and M.D. 
Lukin, Phys. Rev. Lett. 89, 220407 (2002). 

9 M. Kohl, H. Moritz, T. Stoferle, K. Giinter and T. 
Esslinger, Phys. Rev. Lett. 94, 080403 (2005). 

10 R. Jordens, N. Strohmaier, K. Giinter, H. Moritz and T. 
Esslinger, Nature 455, 204 (2008). 

11 U. Schneider, L. Hacker miiller, S. Will, Th. Best, I. Bloch, 
T.A. Costi, R.W. Helmes, D. Rasch and A. Rosch, Science 
322, 1520 (2008). 

12 J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, P 
Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer and 
A. Aspect, Nature 453, 891 (2008). 

13 M. White, M. Pasienski, D. McKay, S.Q. Zhou, D. Ceper- 
ley and B. DeMarco, Phys. Rev. Lett. 102, 055301 (2009). 

14 L. Fallani, J.E. Lye, V. Guarrera, C. Fort and M. Inguscio, 
Phys. Rev. Lett. 98, 130404 (2007). 

15 K. Giinter, T. Stoferle, H. Moritz, M. Kohl and T. 
Esslinger, Phys. Rev. Lett. 96, 180402 (2006). 

16 S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, 
K. Sengstock and K. Bongs, Phys. Rev. Lett. 96, 180403 
(2006) 

17 M. Theis, G. Thalhammer, K. Winkler, M. Hellwig, G. 
Ruff, R. Grimm and J. H. Denschlag, Phys. Rev. Lett. 93, 
123001 (2004). 

18 J.T. Stewart, J.P Gaebler and D.S. Jin, Nature 454, 744 
(2008). 

19 Q. Chen and K. Levin, Phys. Rev. Lett. 102, 190402 



(2009). 

20 A. Damascelli, Z.-X. Shen, Z. Hussain, Rev. Mod. Phys. 
75, 473 (2003). 

21 W. Metzner and D. Vollhardt, Phys. Rev. Lett. 62, 324 
(1989). 

22 D. Vollhardt, Correlated Electron Systems, vol. 9, ed. V. J. 
Emery, (World Scientific, Singapore, 1993), p. 57. 

23 A. Georges, G. Kotliar, W. Krauth and M.J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

24 G. Kotliar, S.Y. Savrasov, K. Haule, V.S. Oudovenko, O. 
Parcolletm and C.A. Marianetti, Rev. Mod. Phys. 78, 865 
(2006). 

25 R. Vlaming and D. Vollhardt, Phys. Rev. B 45, 4637 
(1992). 

26 M. Ulmke, V. Janis and D. Vollhardt, Phys. Rev. B 51, 
10411 (1995). 

27 V. Dobrosavljevic and G. Kotliar, Phys. Rev. Lett. 78, 
3943 (1997); Phil. Trans. R. Soc. Lond. A 356, 57 (1998). 

28 Y. Song, R. Wortis and W.A. Atkinson, Phys. Rev. B 77, 
054202 (2008). 

29 V. Dobrosavljevic, A. A. Pastor and B.K. Nikolic, Euro- 
phys. Lett. 62, 76 (2003). 

30 K. Byczuk, W. Hofstetter and D. Vollhardt, Phys. Rev. 
Lett. 94, 056404 (2005). 

31 K. Byczuk, Phys. Rev. B 71, 205105 (2005). 

32 M.C.O. Aguiar, V. Dobrosavljevic, E. Abrahams, and G. 
Kotliar, Phys. Rev. Lett. 102, 156402 (2009). 

33 K. Byczuk, W. Hofstetter and D. Vollhardt, Phys. Rev. 
Lett. 102, 146403 (2009). 

34 E. Miranda and V. Dobrosavljevic, Phys. Rev. Lett. 86, 
264 (2001). 

35 E. Miranda and V. Dobrosavljevic, J. Magn. Magn. Mater. 
226-230, 110 (2001). 

36 M.C.O. Aguiar, E. Miranda and V. Dobrosavljevic, Phys. 
Rev. B 68, 125104 (2003). 

37 R. Abou-Chacra, P.W. Anderson and D.J. Thouless, J. 
Phys. C 6, 1734 (1973). 

38 A. Alvermann and H. Fehske, J. Phys.: Conf. Ser. 35, 145 
(2006). 

39 K. Byczuk, W. Hofstetter and D. Vollhardt, Phys. Rev. B 



10 



69, 045112 (2004). 

40 M. Potthoff, T. Wegner and W. Nolting, Phys. Rev. B 55, 
16132 (1997). 

41 H. Kajueter and G. Kotliar, Phys. Rev. Lett 77, 131 
(1996). 

42 G. Kotliar and A.E. Ruckenstein, Phys. Rev. Lett. 57, 1362 
(1986). 

43 X.Y. Zhang, M.J. Rozenberg and G. Kotliar, Phys. Rev. 
Lett. 70, 1666 (1993). 

44 U. Yu, K. Byczuk and D. Vollhardt, Phys. Rev. Lett. 100, 
246401 (2008). 

45 S. Kirkpatrick, B. Velicky and H. Ehrenreich, Phys. Rev. 
B 1, 3250 (1970). 

46 A. Gonis, Green's functions for ordered and disordered sys- 
tems, (North- Holland 1992). 

47 K. Byczuk, M. Ulmke and D. Vollhardt Phys. Rev. Lett. 
90, 196403 (2003). 

48 If the random position of the immobile atoms are sub- 
ject to thermal fluctuations, so called annealed disorder, 
the proper microscopic Hamiltonian should be the Falicov- 
Kimbal model, where mobile and immobile components are 
thermodynamically coupled via on-site interaction, cf. J.K. 
Freericks and V. Zlatic, Rev. Mod. Phys. 75, 1333 (2003). 

49 E.N. Economou, Greens Functions in Quantum Physics 
(Springer- Verlag, Berlin, 3rd ed., 2006), p. 357. 

50 M. Eckstein, M. Kollar, K. Byczuk, and D. Vollhardt, 
Phys. Rev. B 71, 235119 (2005). 

51 M. Kollar, M. Eckstein, K. Byczuk, N. Blumer, P. van Don- 
gen, M.H. Radke de Cuba, W. Metzner, D. Tanaskovic, V. 
Dobrosavljevic, G. Kotliar, and D. Vollhardt, Ann. Phys. 
(Leipzig) 14, 642 (2005). 

52 In order to minimize the statistical fluctations we artiflcally 



increase the size of our ensemble after having reached the 
self-consistency. This is done by constructing the histogram 
on basis of typically 50—100 successional update iterations. 

53 E. Miiller-Hartmann, Z. Phys. B 76, 211 (1989). 

54 P. Soven, Phys. Rev. 156, 809 (1967). 

55 R.J. Elliott, J. A. Krumhansl and P. L. Leath, Rev. Mod. 
Phys. 46, 465 (1974). 

56 D. Lloyd, J. Phys. C 2, 1717 (1969). 

57 D. Thouless, Phys. Reports 13, 93 (1974). 

58 F. Wegner, Z. Phys. B 44, 9 (1981). 

59 G.R. Reich and P.L. Leath, J. Stat. Phys. 19, 611 (1978). 

60 E.N. Economou and M.H. Cohen, Phys. Rev. B 5, 2931 
(1972). 

61 S. Kirkpatrick and T.P. Eggarter, Phys. Rev. B 6, 3598 
(1972). 

62 A. Alvermann and H. Fehske, Eur. Phys. J. B 48, 295 
(2005). 

63 G. Schubert, A. WeiBe and H. Fehske, Phys. Rev. B 71, 
045126 (2005). 

64 P. Lombardo, R. Hayn and G.I. Japaridze, Phys. Rev. B 
74, 085116 (2006). 

65 D. Semmler, PhD thesis, Frankfurt University, 2010. 

66 J.F.L. Hopkinson, J.B. Pendry and D.J. Titterington, 
Comput. Phys. Commun. 19, 69 (1980). 

67 P.J. Durham, J. Phys. F 11, 2475 (1981). 

68 J. Kudrnovsky and J. Masek, J. Phys. F 16, 943 (1986). 

69 A. Bansil and M. Pessa, Phys. Scr. T4, 52 (1983) and 
references therein. 

70 Y. He , Q.J. Chen and K. Levin, Phys. Rev. A 72, 
011602(R) (2005). 



